Thermodynamics of doped Kondo insulator in one-dimension 
- Finite-temperature DMRG study - 
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The finite-temperature density-matrix renormalization-group method is applied to the one- 
dimensional Kondo lattice model near half filling to study its thermodynamics. The spin and charge 
susceptibilities and entropy are calculated down to T = O.OSt. We find two crossover temperatures 
near half filling. The higher crossover temperature continuously connects to the spin gap at half 
filling, and the susceptibilities are suppressed around this temperature. At low temperatures, the 
susceptibilities increase again with decreasing temperature when doping is finite. We confirm that 
they finally approach to the values obtained in the Tomonaga-Luttinger (TL) liquid ground state for 
several parameters. The crossover temperature to the TL liquid is a new energy scale determined 
by gapless excitations of the TL liquid. The transition from the metallic phase to the insulating 
phase is accompanied by the vanishing of the lower crossover temperature. 



The Kondo insulator is a typical strongly correlated in- 
sulator and develops a spin gap at low temperatures [Q . 
The half-filled Kondo lattice (KL) model has been stud- 
ied as its theoretical model, particularly in one dimen- 
sion (ID) by both numerical and analytical approaches 
. The ground state of this model is shown to have 
both spin and charge gaps, Ag, Ac, and strong correla- 
tion effects appear in their difference Q. The spin gap 
is always smaller than the charge gap, and for small ex- 
change coupling J the spin gap is exponentially small. 
As ~ 6xp( -l/apJ) 1,1^, while the charge gap is lin- 
ear in J 1,1. 

At finite temperatures, the spin gap characterizes 
unique temperature dependence of the excitation spec- 
trum. With increasing temperature, the structure of the 
charge gap in the dynamic charge structure factor and 
the quasiparticle density of states disappears at T As, 
much lower than Ac, although the latter is the energy 
scale of charge excitations at T = [^|j9j. This feature 
is also seen in the temperature dependence of the charge 
susceptibility, which drastically decreases below T ~ Ag 
pof . As for the spin susceptibility, it decreases exponen- 
tially with the energy scale of Ag , as expected . 

When a finite density of carriers are doped, the ID KL 
model belongs to another universality class. The ground 
state and the low energy excitations are described as a 
Tomonaga-Luttinger (TL) hquid In contrast to 

the half filled case, the ground state has gapless excita- 
tions in both spin and charge channels. Consequently 
the spin and charge susceptibilities are finite at T = 
and determined by the velocities of the collective excita- 
tions, Xs = l/(27r?;s) and Xc = '2Kc/ {ttVc) where is 
the Luttinger- liquid parameter whereas Xs = 0, 

Xc = at half filling due to the gap. 

In the present paper we study thermodynamics of the 
ID KL model near half filling, i.e. in the vicinity of 
metal-insulator transition. Thermodynamic quantities 
are calculated by the finite-temperature density-matrix 
renormalization-group (finite-T DMRG) method [p2|-p5[ , 
and we find that the spin gap at half filling appears as a 



crossover temperature such that the susceptibilities are 
suppressed around there. We also find an new lower 
crossover temperature which is the energy scale of the 
gapless excitations of the TL liquid ground state. This is 
sensitive to the hole doping, and the transition from the 
metallic phase to the insulating phase corresponds to the 
vanishing of the lower crossover temperature. 

The Hamiltonian we use in the present study is the ID 
KL model described as 

T-i = -<^(4^Ci+is + H.C.) + J ^ S,; • ^CTss'cl^Cis', 

i,s i,s,s' 

(1) 

where cTss' are the Pauli matrices and Si = 
s' l'^ ss' fisfis' is the localized spin at site i. The 
model has hoppings — t [t > 0) for only nearest neigh- 
bor pairs. The density of conduction electron ric is unity 
at half filling, and hole doping (ric = 1 — 5) is physically 
equivalent to electron doping (tic = 1 + 6) due to the 
particle-hole symmetry. 

In order to study thermodynamics we employ the 
finite-T DMRG method (2|-|2|]. By iteratively increas- 
ing the Trotter number of the quantum transfer matrix, 
we can obtain the eigenvector for the largest eigenvalue 
with desired accuracy. Thermodynamic quantities are 
directly calculated from this eigenvector, and the extrap- 
olation in the system size is not needed . This method 
was first applied to the quantum spin systems and shown 
to be reliable down to the low temperature T = O.OIJ 
p^ , ^ . This method is free from statistical errors and 
the negative sign problem, which are advantages com- 
pared with the quantum Monte Carlo method. 

To obtain thermodynamic quantities with fixed hole 
density (5 = 1 — ric , we need chemical potential ^ at each 
temperature. This requires many DMRG calculations 
at different fixed chemical potentials: for each J we use 
36 sets of n with a typical interval = 0.025t. The 
number of states used in the present study is typically 54 
and corresponding truncation error is 10"'^ at the lowest 
temperature T = O.OM with the Trotter number 60. 
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FIG. 1. Chemical potential of the one-dimensional Kondo 
lattice model near half-filling. 

The T-dependence of the chemical potential for several 
(5's is shown in Fig. 1 for J = 1.6t and 1.2t. Note that 
the chemical potential at half filling is always zero due 
to the particle-hole symmetry. At high temperatures, fi 
is similar to the value for the free conduction electrons 
and indicates metallic behavior. At low temperatures, 
however, a significant increase in | — /i| appears for small 6. 
In the limit of T = 0, /i approaches T = value and this 
clearly shows the presence of the quasiparticle gap Aqp 
as shown in Fig. 1. Previous T — DMRG calculations 
show that Aqp = 0.7t for J 1.6t and 0A7t for J = 
1.2i, which is consistent with the present calculation. In 
the following, we calculate thermodynamic quantities for 
fixed /i's and convert them for fixed (5's using these data 
ofM'5,T). 

Temperature dependence of spin susceptibility Xs{T) 
is plotted in Fig. 2 for J/t = 1.6 and 1.2 at < (5 < 0.2. 
At high temperatures Xs is asymptotically determined by 
the sum of the Curie term of the localized spins, l/(4r), 
and the Pauli susceptibility, XPauU) of the free conduction 
electrons. This is actually seen in the inset of Fig. 2(a). 
For 5 < 0.2, the change in the density of states of the free 
conduction band due to hole doping is within 5% leading 
to little change in XPauU, and it is reasonable that the 
(5-dependence of Xs is small at high temperatures. 

With decreasing temperature the spin susceptibility in- 
creases owing to the Curie term of the localized spins, 
but around T ~ Ag, the spin susceptibility starts to 
be suppressed as in the (5 = case. Previous T — 
DMRG calculations show that A^ = QAt for J = IM, 
and As = 0A6t for J — l.2t These behaviors suggest 
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FIG. 2. Spin susceptibility of the one-dimensional Kondo 
lattice model near half-filling. We use the units of (gfJ-B)^ ~ 1- 



that the spin gap at S — persists as the crossover tem- 
perature characterizing the suppression of Xs even away 
from half filling. 

With further decreasing temperature, Xs sharply in- 
creases again when doping is finite, whereas it exponen- 
tially decreases at half filling with the energy scale of Ag. 
The increase in Xs seems to be proportional to S at low 
temperatures. In order to see this (5-dependence in more 
detail, we plot the difference in Xs between S > and 
(5 = divided by S. As shown in Fig. 3 for J = 1.6t, 
the universal behavior is observed at low temperatures 
indicating Xs(T) - Xs{T,6 = 0) + (5/(4T). This means 
that the doped holes induce almost free spins of 5 = i 
with density S. 

In the limit of T — > 0, the thermodynamic properties 
are determined by the gapless collective excitations of the 
TL liquid, and the spin susceptibility is given by the spin 
velocity Vs as Xs = l/(27r?;s) ||l^ , M . Thus there must 
be a crossover temperature where Xs{T) deviates from 
Xs{T) - Xs(T,(5 = 0) +(5/(4r) towards l/(27r?;s). This 
crossover temperature is expected to be proportional to 
Vs, because it is the energy scale of the spin excitations. 
One can estimate this crossover temperature from Xs at 
T = through the relation Vg = l/(27rxs)- 

The Xs(r=0) are plotted in Fig. 4 for J = 2.0t. They 
are calculated by the zero-temperature DMRG method 
|p2| with open boundary conditions via the size depen- 
dence of the lowest spin excitation energy, AE{L) = 
TTVs/L. The T = susceptibility is very sensitive to the 
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FIG. 3. Low temperature behavior of spin susceptibility 
near half-filling. 
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FIG. 4. Doping dependence of spin susceptibility at T = 0. 

hole doping 6 and diverges in the Umit of 5 — > 0. It 
reflects the behavior that the spin velocity as a charac- 
teristic energy scale vanishes. Thus the crossover temper- 
ature to the TL liquid phase may correspondingly vanish 
as (5 ^ 0. The ^-dependence of Xa seems to be exponen- 
tial in the present calculation, but we will discuss this 
point again later. 

The crossover to the TL liquid is more clearly seen in 
the charge susceptibility, Xc- The charge susceptibility is 
defined as the change in the conduction electron density 
due to a small shift of chemical potential, Xc = duc/dii — 
—dS/d^. The results for J = l.Qt and 1.2t are shown in 
Fig. 5. 

At high temperatures, Xc is almost the same as that of 
free conduction electrons and proportional to 1/T. With 
decreasing temperature, Xc for small doping is suppressed 
around T ~ Ag, as for Xs- The suppression of Xc at 
(5 = is due to the development of the charge gap. Below 
T ~ As, the correlation length of the localized spins be- 
comes longer than the charge correlation length, and this 
induces an internal staggered magnetic field for the con- 
duction electrons through the Kondo coupling. When the 
charge excitations are concerned, their time scale is much 
shorter than that for the spin excitations (oc A~^), and 
the staggered field may be considered as almost static. 
This staggered field induces a unit-cell doubling for the 
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FIG. 5. Charge susceptibility of the one-dimensional 
Kondo lattice model near half-filling. The dot on the vertical 
axis is the T — value calculated by the zero-temperature 
DMRG method. 

conduction electrons and a finite gap appears at the 
Fermi energy of the otherwise free conduction electrons. 
Since this behavior is induced by the development of long- 
range spin correlations, it becomes noticeable only below 
T ~ As rather than Ac. The present results show that 
contributions to Xc from the small number of holes are 
small at temperature T As and the suppression is vis- 
ible near half filling. 

With further decreasing temperature, the temperature 
dependence of Xc at (5 > differs from the 6 = case 
as for Xs{T). Xc at half filhng becomes exponentially 
small below T ~ As with the energy scale of the quasi- 
particle gap, Xc ~ exp(— Aqp/T) but for finite 6, Xc 
sharply increases. The increase of Xc at low temperatures 
seems to be larger for smaller S. This behavior near half 
filling is consistent with the doping dependence of Xc at 
T = 0, which is obtained by the zero-temperature DMRG 
method. At T = 0, Xc is given by the difference in the 
chemical potential Xc = L{n{L, Nc + 1) — fi{L, Nc — l)), 
where 2//(L,iVc + 1) = Eg{L,Nc + 2) - Eg{L,Nc) and 
Eg{L, N) is the ground state energy of N conduction elec- 
trons in the system of length L. The results are shown in 
Fig. 6, which actually shows that the Xc increases with 
decreasing 6 and it seems to diverge in the limit of (5 — > 0. 
The divergent behavior of Xc corresponds to the vanish- 
ing charge velocity as the characteristic energy of the TL 
liquid as (5 — > 0. The i5-dependence of Xc is close to l/S 
and different from the exponential dependence observed 
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spin gap at half filling and the other is the characteristic 
energy of collective excitations in the TL-liquid ground 
state. The energy scales of the TL liquid becomes smaller 
with approaching half filling and seems to vanish in the 
limit of (5 — + 0. 
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from the Ministry of Education, Science, Sports and Cul- 
ture of Japan. 
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FIG. 6. Doping dependence of charge susceptibility at 
T = 0. The broken line is a guide for eyes. 
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FIG. 7. Entropy of the one-dimensional Kondo lattice 
model near half filling. 

in Xs in the present DMRG calculation. The Xs seems 
to diverge stronger than Xc but this point is not clear at 
the moment. 

The new energy scale characterizing gapless excitations 
of the TL liquid may also be seen in the temperature de- 
pendence of entropy S(T) . The entropy is obtained as the 
derivative of the free energy with respect to temperature, 
and the results are shown in Fig. 7 for J = 1.6t. At low 
temperatures the entropy at half filling becomes exponen- 
tially small reflecting finite spin and charge gaps. When 
6 > 0, the entropy is enhanced at low temperatures. This 
enhancement is due to the appearance of gapless excita- 
tions in the TL liquid. Entropy of the TL liquid is pro- 
portional to T at low temperatures, S = TrT{v~^ +v~^)/3 
[po| , pT| . However, for small 6, the T-dependence is small 
even at temperatures T ~ O.OSt. The remaining entropy 
is consistent with the entropy of free spin-i carriers with 
density 6. This implies that the characteristic energy 
scale for gapless excitations of the TL liquid is small near 
half filling and the T-linear dependence will be seen at 
further low temperatures T ^ O.OSt. 

In conclusion, we have applied finite-T DMRG method 
to the one-dimensional Kondo lattice model away from 
half filling, and found two crossover temperatures which 
characterize thermodynamics near half filling. One is the 
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